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ABSTRACT 


The  FPS-16  tracking  systems  superimpose  undesired  oscillations 
on  the  real  position  data  resulting  in  rapidly  fluctuating  successiv 
position  points  which  are  physically  unrealistic.  A  linear  digital 

M*S8 

filter  of  the  form  (L  «  £  W  Q  is  applied  separately  to  each 

*  M--58 

component  to  smooth  the  data.  The  frequency  response  is  given  and  ? 
data  are  corrected  by  a  method  derived  by  Eddy  et  al.  (196S),  Un¬ 
desired  high  frequency  oscillations  are  effectively  eliminated  and 
successive  wind  profiles  show  good  continuity. 
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INTRODUCTION 


Atmospheric  temperature,  pressure,  density,  and  wind  data, 
derived  from  meteorological  rocket  soundings  are  used  in  meteor¬ 
ological  research  and  for  military  projects  requiring  a  knowledge 
of  stratospheric  behavior.  Detailed  wind  profiles  are  often  required 
for  computing  trajectories  for  ballistic  rockets.  This  paper  will 
describe  the  filtering  and  correction  techniques  used  to  derive  winds 
from  an  FPS-16  radar  track  of  a  falling  object. 


FILTERING  TECHNIQUE 


Wind  data  derived  from  radar  tracks  of  parachutes  and  spheres 
falling  through  the  atmosphere  are  used  as  meteorological  support 
data  by  several  projects  at  White  Sands  Missile  Range,  New  Mexico. 

An  FPS-16  radar  tracks  the  sensor  and  records  its  position 
relative  to  the  radar  on  a  magnetic  tape  at  the  rate  of  twenty 
points  per  second.  The  meteorological  wind  reduction  technique, 
however,  uses  only  every  other  point,  or  ten  points  per  second. 

These  position  points  are  specified  by  time  and  three  space  coot 
dinates  (slant  range,  azimuth  angle,  and  elevation  angle).  A  correc¬ 
tion  for  refraction  and  earth  curvature  is  then  applied  to  the 
position  points.  A  typical  wind  profile  consists  of  about  18.00C 
data  points  and  yields  wind  data  from  25  km  to  about  65  km. 

As  is  true  with  any  tracking  system,  the  system  itself  super¬ 
imposes  undesired  oscillations  (noise)  on  the  real  position  data 
resulting  in  rapidly  fluctuating  successive  position  points  which 
are  physically  unrealistic.  To  help  compensate  for  this  feature 
a  linear  digital  filter,  eq  (1),  in  the  form  of  a  weighted  running 
average  over  117  points  *(11.7  seconds)  is  applied  separately  to 
each  component  (slant  range,  azimuth  angle,  and  elevation  angle). 

M-S8 

^  "  M.-58  KmQ|C*M.  (1) 

where  WM  is  the  value  of  the  Mth  weight,  is  the  (K«-M)th  un¬ 

smoothed  value  of  a  coordinate,  and  is  the  Kth  smoothed  value 

of  the  coordinate.  The  weights  are  symmetrically  centered  about 
Wo^WN  *  W~N^  *nd  ar®  shown  in  Figure  1.  Assuming  the  unsmoothed 

coordinate  values  may  be  represented  by  a  sum  of  sinusoidal  os¬ 
cillations  of  various  amplitudes,  phases,  and  frequencies,  a  fre¬ 
quency  response  (ratio  of  the  amplitude  of  a  sinusoidal  wave  in 

*  Table  1  shows  the  coefficients  for  the  117-point  symmetrical  filter. 
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TABLE  1 


WGT(O) 
WGT(l) 
WGT  (2) 
WGT(3) 
WGT (4) 
WGT (5) 
WGT (6) 
K'CT(V) 
WGT (8) 
WGT (9) 
WGT (10) 
WGT (11) 
WGT(1 2) 
WGT (13) 
WGT (14) 
WGT(15) 
WGT (16) 
WGT (17) 
WGT (IS) 
WGT (19) 
WGT (20) 
WGT (21) 
WGT (2 2) 
KGT(23) 
WGT(24) 
WGT(2S) 
WGT  (26). 
WGT  (27). 
WCT(28). 
KGT(29). 


COEFFICIENTS  FOR  THE  117-POINr  SYMMETRICAL  FILTER. 


0.3S708064E-01 

WGT (30) 

at 

-0. 10131001E-02 

0.35626344E-01 

WGT(31) 

a 

-0.14721988E-02 

0.3527267SE-01 

WGT  (3  2) 

■ 

-0. 191 44 034 E- 02 

0.34784139E-01 

WGT(33) 

n 

-0.21881496E-02 

0.34025734E-01 

WGT (34) 

M 

-0. 24564 068E-07 

0.33167716E-01 

WGT (35) 

m 

-0.2S690121E-02 

0.3206788SE-01 

KGT(36) 

m 

-0. 26900660E-02 

0.30898912E-01 

WGT(37) 

M 

-0. 26702617E-02 

0.29S20S9SE-01 

WGT (38) 

m 

-0. 26742817E-02 

0.28107178E-01 

WGT(39) 

m 

-0.25S31311E-02 

0.26S19604E-01 

WGT (40) 

m 

-0. 247154 28E-02 

0. 24932820E-01 

WGT (41) 

m 

-0. 22803228E-02 

0. 23 208 04 SE- 01 

WGT (4 2) 

m 

-0. 2143771 0E-02 

0. 21S20080E-01 

WGT (43) 

m 

-0. 19120444E-02 

0. 19729S87E-01 

WGT(44) 

m 

-0. 17486000E-02 

0. 18010426E-01 

WGT  (45) 

m 

-0. 15025814E-02 

0, 16221950E-01 

WGT (46) 

at 

-0. 13362864E-02 

0.14S36376E-01 

WGT(47) 

m 

-0.10975813E-02 

0.1 281 1 107E-01 

WGT(48) 

a 

-0.94739597E-03 

0.1 1 21614SE-01 

WGT (49) 

«j 

-0.7321 66 18E- 03 

0.96064442E-02 

WGT (50) 

m 

-0.61135714E-03 

0.81493637E-02 

WGT (51) 

m 

-0. 42993361 E-03 

0.66970661E-02 

WGT(S2) 

m 

- 0. 34 591 034 E- 03 

0.54140097E-02 

WGT(53) 

m 

-  0. 20284763E-03 

0.41493890E-02 

WCT(54) 

m 

-  0. 15742585E-03 

0.3 064 6684 E-0? 

WGT(55) 

a 

- 0.51962753E-04 

0. 20060650E-02 

WGT (56) 

8 

- 0.42007490E-04 

0.1132119SE-02 

WGT  (5  7) 

- 

0. 31 125163E-04 

0.28621851E-03 

-0.37580166E-03 

WGT (58) 

= 

0. 12746845E-04 

f 


WEIGHT  NUMBER  - 


the  smoothed  data  to  the  amplitude  of  the  same  wave  in  the  unsmoothed 
data)  may  be  calculated,  eq  (2).  This  will  give  the  effect  of  the 
filter  on  the  data. 

M-S8  . 

R(f)  -  l  wu  COS  (iiil)  for  0  <  f  <  S  sec-1,  (2) 
M—  58  M  * 

where  R(f)  is  the  frequency  response  at  frequency1  f.  The  frequency 
response  of  this  filter  is  shown  in  Figure  2.  The  filter  essentially 
eliminates  oscillations  which  have  frequencies  greater  than  0.3  sec-' 
(periods  less  than  three  seconds).  Figures  3,  4,  and  S  show  typical 
samples  of  the  first  differences  in  the  smoothed  and  unsmoothed 
values  of  slant  range,  azimuth,  and  elevation  angles  at  increments 
of  0.1  second. 


The  first  and  second  derivatives  of  each  coordinate  are  approxi¬ 
mated  by  equations  (3)  and  (4). 


*  ^K-p1 


ioo(?L 


-  2\  *  \ 


where  QR,  and  QR  are  the  Kth  values  of  the  smoothed  coordinate 

and  its  first  and  second  derivatives,  respectively.  A  transforma¬ 
tion  of  coordinates  is  performed  which  gives  position,  velocity, 
and  acceleration  data  in  terms  of  components  oriented  north-south 
(y),  east-west  (x),  and  normal  to  the  surface  of  the  earth  (z). 

The  acceleration  values  of  each  component  fluctuate  excessively 
and  are  physically  unacceptable.  Ten  weights  are  used  to  filter 
these  data  and  are  derived  from  equation  (S), 

\  *  m£0  V  K-10M,  (S 

where  A^  is  the  Kth  value  of  the  smoothed  acceleration  data, 
the  Mth  weight,  and  A^,  the  (K-lOM)th  value  of  the  unsmoothed 

acceleration  data  of  one  of  the  components  .  Unsmoothed  values 
one  second  apart  instead  of  a  tenth  of  a  second  apart  are  used  in 
the  smoothing.  The  amplitude  of  the  frequency  response  may  be 
calculated  from  equations  (6),  (7),  and  (8). 


|R*(f)|  -  (C[f]2  ♦  S(fJ2)1/2, 


C(f)  -  1  cos  (2wfM), 

M-0 


*  Table  2  shows  the  coefficients  for  the  10-point  non-symmetrical 
filter. 
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COEFFICIENTS  FOR  THE  10- POINT  NON -SYMMETRICAL  FILTEP 

KGT(O)  -  0.93263928E-01 
WCT(I)  *■  0. 10963999E  00 
’,.'GT(2)  »  0.1214S806E  00 
fc'GT (3)  =  0.1278C346E  00 
WGT (4)  =  0.12803132E  00 
h'GT(S)  =  0. 12177899E  00 
WGT (6)  «  0. 10896906E  00 
KGT(7)  -  0. 89803061E-01 
KGT(8]I  =*  0.64 7466S8E-01 
KGT(9)  -  0. 34506476E-01 


S(f)  -  l  K  sin  (2*fM) 

-1  M“°  • 

for  0  <  f  <  0.5  sec  where  |R  (f) |  is  the  magnitude  of  the  fre- 
quency- response  at  frequency  f.  The  weights  are  shown  in  Figure 
6  and  |R*(f)|  in  Figure  7. 


(8) 


CORRECTION  TECHNIQUE 


Another  error  of  the  system  must  be  corrected,  namely  the 
sensor's  ability  to  respond  to  the  actual  wind.  The  faster  the 
sensor  falls,  the  less  likely  it  will  respond  to  the  actual  wind. 

In  other  words,  small-scale  wind  oscillations  will  have  little 
effect  on  the  sensor  while  those  with  longer  periods  will  be  ob¬ 
served  with  greater  accuracy.  Therefore,  the  wind  sensor  itself 
becomes  a  time  varying  filter  applied  to  the  wind  data.  Eddy  ct  al. 
(1965)  designed  a  correction  technique  which  theoretically  eliminates 
this  effect,  eq  (9) . 


V  XK  -  - 


XKZK 


zK  ♦  g 


(9) 


where  Tj,.  X X^,  ZK,  and  ZR  are  the  Kth  values  of  the  east-west 


wind  component,  sensor  velocity,  sensor  acceleration,  sensor  vertical 
velocity,  and  sensor  vertical  acceleration,  respectively,  and  g  the 
gravity  constant.  An  analogous  equation  is  also  applied  to  the 

north-south  component.  is  always  assumed  to  be  zero.  This 

correction  technique  has  been  experimentally  verified  by  Kays  and 
Olsen  (1966).  Typical  north-south  and  east-west  wind  component 
profiles  in  final  filtered  and  corrected  form  are  shown  in  Figure  8. 


CONCLUSION 


The  wind  reduction  technique  discussed  is  most  applicable  when  the 
sensor  is  a  point  source  and  is  tracked  by  a  radar  capable  of  at 
least  a  ten  points  per  second  sampling  rate.  Undesired  high  frequency 
oscillations  are  effectively  eliminated,  and  successive  wind  profiles 
show  good  continuity.  This  technique  is  currently  used  at  White  Sands 
Missile  Range,  New  Mexico,  by  the  Atmospheric  Sciences  Laboratory. 
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The  FPS-16  tracking  systems  superimpose  undesired  oscillations  on  the  real 
position  data  resulting  in  rapidly  fluctuating  successive  position  points  which 
are  physically  unrealistic.  A  linear  digital  filter  of  the  form 
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is  applied  separately  to  each  component  to  smooth  the  data.  The  frequency  response 
is  given  and  the  data  are  corrected  by  a  method  derived  by  Lddy  ct  al.  ( 1965t  • 
Undesired  high  frequency  oscillations  are  effectively  eliminated  and  successive 
wind  profiles  show  good  continuity. 
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